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Abstract 

The quantum elliptic 77-matrices of Baxter-Belavin type satisfy the associative Yang- 
Baxter equation in Mat(IV, C)® 3 . The latter can be considered as noncommutative ana¬ 
logue of the Fay identity for the scalar Kronecker function. In this paper we extend the list 
of 77-matrix valued analogues of elliptic function identities. In particular, we propose coun¬ 
terparts of the Fay identities in Mat (A, C)® 2 . As an application we construct 77-matrix 
valued 2 N 2 X 2 N 2 Lax pairs for the Painleve VI equation (in elliptic form) with four free 
constants using Z^r x Z^v elliptic 77-matrix. More precisely, the four free constants case 
appears for an odd N while even TV’s correspond to a single constant. 
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1 Introduction and summary 


In this paper we continue the study of identities for quantum (and classical) /?-matrices, which 
are similar to the elliptic functions identities for scalar elliptic functions [T3J I5j- More concretely, 
we prove the Fay identities in Mat (N, C)® 2 . It allows us to construct multidimensional Lax pairs 
for the Painleve VI equation with the /?-matrices as matrix elements. 

We start with the list of properties and identities for elliptic functions, and then give their 
/?-matrix version. Most of the properties are known from [2111], HU, [3UI5], [13] and |S]. 

Consider the following functions: 


(p(z,u) = 


d'(0 + u ) 


Ei(z) = 


,V(z) 
$(z) 


1 '{)"'({)) 

E 2 (z) = -d gEl {z) = p(z) - , 


where $(z) is the odd Riemann theta-function 

$(z) = $(z\t) = exp (mT(k + -) 2 + 27r*(2: -f— ){k -|—) j 
kez ^ ' 


( 1 . 1 ) 

( 1 . 2 ) 


(1.3) 


and p(z) is the Weierstrass ^-function. 

Following [16] the function (11.11) is referred to as the Kronecker function, and (II.2ft are called 
the (first and the second) Eisenstein functions. 

The Kronecker function can be considered as a section of the Poincare bundle V over E r x E(_. 
Here E T is the elliptic curve 


E t = C/(Z + rZ), 5 suit > 0 , 


(1.4) 


is its Jacobian 


(E ' T ~ E t ). The Poincare bundle V is a line bundle over E r x E' r 
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V 


( 1 . 5 ) 


S T x S' 



£ r K 

specialized by (11.61) . (1 1.71) . (ll.lOp and (11.lip . 

The properties of theta-function (jl.3p (including Riemann identities, see jTT] ) provides the 
following set of properties and relations for the functions (U.1[) - (I1.2I) : 

• Arguments symmetry: 


4>(z, u ) = 4>(u , z) , z G S T , u G S(_, 


( 1 . 6 ) 


• Local expansion: 

4(z,u) = 1 + Bi(«) + Z -(El(n) - p(n)) + 0(z 2 ), (1.7) 


• Residues: 


Res (j){z , u ) = Res 4>(z,u) = Res E x (z) = 1, 

z=0 u =0 2=0 


• Parity: 

(j)(-z, -u) = -<j>(z,u) , #i(-,z) = --Eu(z), E 2 (-z) = E 2 (z) , 

• (Quasi)periodicity properties: 

4>(z + l,u) = (p(z,u), Ei(z + 1) = Ei(z), E 2 (z + 1) = E 2 (z) , 

4>(z + r, u) = e~ 2mu (j)(z, u ), Td ( 2 ; + r) = Td (z) - 2m , E 2 (z + r) = E 2 (z ), 

• i/eat equation: 

27 nd T (j)(z, u ) = d z d u (j)(z, u ), 


• Derivatives: 

d u </>(z,u) = 4>(z,u)(Ei(z + u) - E^u )), 
d z 4>(z,u) = <t>(z,u)(E 1 (z + u) - E^z)), 

• Fay (trisecant) identity JE|/; 

<^>(a, u)<f>(y, w) — <^>(a — y, a)0(y, a + w) + )>(y — x, w)<f>(x, u + ui), 


( 1 . 8 ) 

(1.9) 


( 1 . 10 ) 

( 1 . 11 ) 

( 1 . 12 ) 

(1.13) 

(1.14) 


(1.15) 
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Degenerated Fay identities: 


<f>(x, z)<f{x,w) = (j)(x, z + w)(E 1 (x) + E 1 (z) + E 1 (w) - E^x + z + w)), (1.16) 

or 

z)(f>(y, z) = (j>(x + y, z)(E i(a?) + E x (y) + E^z) - E x {x + y + z )), (1.17) 


4>(x, z)(fi(x, -z) = E 2 (x) - E 2 (z) = p(x) - p(z) . (1.18) 

• Geometric interpretation: The Kronecker function (f)(z, u ) is a section of the Poincare 
bundle V. It is a line bundle over E T x S r , defined by the conditions (11.6ft . (11.7ft . (11.10ft . 

(imj. 

• Green function: The Kronecker function is the Green function for the operator 8 in the 
space of one forms ^4^ 1 , 0 - ) (S r ) with the boundary conditions (II. 10 ft and (II.lift : 

8(j)(z, u) = 5 2 (z, z) . (1.19) 


Quantum i?-matrices. Consider Zjv x Z^r (Baxter-Belavin’s) elliptic i?-matrix [2J0] in the 
fundamental representation (see also ro- It is defined via the finite-dimensional representation 
of the Heisenberg group: 


27TZ 

Q, A e Mat (IV, C) : Q ki = S ki ex p(—k), A H = < 5 fc -/+i= 0 m Od iv , k,l = l,...,N, 

exp(27Tz^^-)(5 71 A 72 = A 72 Q 71 , 71,72 6 Z . 

Introduce the siu-algebra basis in Mat (AT, C): 


T 7 := T 7l72 = exp(vn^p)Q 71 A 72 , 71,72 = 0,..., N - 1 . 
The same definition is used for any 7 e Z x2 . Then 

T a T p = n a ^T a+p , n a , b = exp (^— (f3 1 oi 2 - p 2 oci)J , 
where a + f3 = (07 + fii,a 2 + ft 2 ). The Akmatrix is defined as 

Ri 2 (u) = co a + h)T a ® T_ q e Mat (AT, C ) 02 , 

wlicrfQ 


TLtjy x it 


<p a (u, co a + h) = exp( 27 nud T ca a )(J)(u, co a + h ), co a = 
The Zat x Z tv symmetry means that for g = Q, A 

0 9 ® g)R i 2 («)(t/ _1 ® 3 _1 ) = ^ 2 («) • 


Oil + Oi 2 T 

N 


1 Here d T oj a = a 2 /N. 


( 1 . 20 ) 

( 1 . 21 ) 

( 1 . 22 ) 

(1.23) 

(1.24) 

(1.25) 

(1.26) 
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For N = 1 the R -matrix (II . 24[) is the scalar Kronecker function cp(h,u) (II. ip . Notice that (jl ,24[) 
is normalized in such a way that the uuitarity condition acquires the form: 

R 12 ( u ) R 2 i(-u) = N 2 (j){Nh, u)(j){Nh, -u) 1 1 = N 2 (p(Nh) - p(u)) 1 ® 1. (1.27) 

The latter can be considered as analogue of 01.181) . Here R 2 i(z) = R i 2 R i 2 ( z ) R i 2 , where 

1 N 

Pl2 = n Ta ® T ~ a = S E v ® E a > ( E a)ki = (1.28) 

a i,j =1 

is the permutation operator. We also use notation Rab( z ) which differs from (I1.27P by T“®T^ a = 
l<g)...l<g)T a <g)l...l<g) T„ a <g) 1... <g) 1 instead of T Q <g> T_ Q (i.e. T a and T__ a are in the a-th and 
6-th components). The number of components in the tensor product is an integer N. It means 
that R^ b is considered as an element of Mat(IV, C) 0 ^, i.e. N^ x N R matrix. 

The properties and identities (ll.8D - fll.17D have the following analogues for i?-matrices: 

• Arguments symmetry: 

R i 2 (z) = R f 2 (Nh)P 12 , 

• Local expansion in h is the classical limit: 

r i 2 (z) = ft~ l 1 ® 1 + 7 - 12 ( 2 ) + hm 12 (z) + 0 (h 2 ), 
where ri 2 {z) is the classical (Belavin-Drinfcld [4]) r-matrix: 

7 - 12 ( 2 ) = Ei(z) 1 <g) 1 + s ^2 / p a (z) T a <g) T_ Q 

Q^O 


(1.29) 

(1.30) 

(1.31) 


and 


77712(2) = 


E 2 (z) - p(z) 


1 + ^ exp(2mzd T uj a )d u (j)(z, u) \ u=u 


T n <g) T_ r 


(1.32) 


a^O 


Similarly to (11.71) we have: 


r, 2 (z) - 2 mi 2 (z) = 1 ® 1 N 2 p(z ), 


(1.33) 


i.e. the quantum /?-matrix is a matrix analogue of the Kronecker function (II.ip while the 
classical one is the analogue of the first Eisensteiu function (1 1.2 j) . 

Expansion with respect to 2 (near 2 = 0) is as follows: 


= — + rT + o(z ), 


wherqj 


rT ] — ^c. ® T-c. (Ei(h + ui a ) + 2 md T u a ) 


(1.34) 


(1.35) 


2 Ri 2 °' > appears as a part of the inverse inertia tensor for relativistic tops [9]. 
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• Residues 


(1.36) 


Res R h V2 (z) — 1 <8> 1 , Res R^Jz) = Res r\ 2 (z) = NP 12 , 
h= o 2=0 2=0 

• Parity: 

R\ 2 {z) = - R 2 i(~ z ) 1 r i20) = -r 2 i(-z ), m 12 (z) = m 21 (-z). (1.37) 

The R-matrix analogue of E 2 {u ) = E 2 (—u ) (H3 appears as E( 2 (^) — — d u r\ 2 (u) (It is 
natural because ri 2 (w) is the analogue of E\{u)). The classical r-matrix is odd. Hence 
F\ 2 {u) is even matrix function. The same answer follows from the local expansions (1 1.71) . 
(| 1 30 | ) : E 2 {u) = - d u <j)(z,u ) | 2=0 , then -d u R{ 2 (u) | 2=0 = -d u r 12 (u). 

• (Quasi)periodicity properties: 


R\ 2 {z + Nut 7 ) = exp(— 2mNhd T uJ 1 ) (T 7 1 <S) 1 )R h l2 (z){T 1 ® 1), 
Rit^iz) = exp(—2mzd r oj J ) ( T~ l <8> l)Rl 2 (z)(l <8) T 7 ). 

In particular, 

R\ 2 {z + 1) = (Q" 1 ® l)R h 12 (z)(Q <8> 1), 

Ri 2 (2 + t) = exp(— 2mh) (A -1 ( 8 ) 1)R^ 2 (^)(A <8) 1), 
R it\z) = R\ 2 (z) , Rit T {z) = exp(—2mz) R r l 2 (z) , 


(1.38) 

(1.39) 

(1.40) 

(1.41) 


r 12 (z + 1) = (Q 1 < 8 > l)ri 2 (^)(Q < 8 > 1), 


r 12 (^ + r) = (A 1 < 8 ) l)r i 2 ( 2 )(A ( 8 ) 1 ) - 2m 1 < 8 > 1 . 
Let us also rewrite (11,39ft as follows: 

R r T /N a. - **) = QpR h M* - *) Q *, 


R h+T/N^ _ ^ = exp(— 27 u ~“ — — ) A„ - 2 t ) As. 

Recall now the R-matrix valued Lax matrix for g $ Calogero-Moser model [ 8 ]: 

N 

C(h) = ^ E a6 (8) Cabih) , C ab {h ) = S ab p a l a <8> l b + 1/(1 - 5 ab )R^ b (z a - Z b ) . 

a,b= 1 


(1.42) 


(1.43) 

(1.44) 


(1.45) 


where E ab is the standard basis of gl^: (E ab ) c d = dac^bd, a,b,c,d = 1...N. Then it follows 
from fjl.43|) - fjl.44jl that 


C(h+1/N) = Q~ 1 C{h)Q, 

C(h + t/N) = exp(-Z/A) \~ l C(R) A exp(Z /N ), 


(1.46) 


where 


N N N 

Q = (J) Qa , A = (J) A a , Z = z a l a 

a= 1 a= 1 a= 1 


(1.47) 


are block diagonal matrices. The number of blocks is N x N, the size of a block is x N- 


tN 
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Heat equation: 


Derivative^: 


2md T R% 2 {z) = d z d h R% 2 (z) 


d h R r l 2 (z) = - (r 12 (z + Nh) Ri 2 (z) + R\ 2 (z) r-&(z - Nh)) 

+Y( E i( z + Nh ) -Mz~ - 2E 1 (Nti)^jR^ 2 (z), 

d z Ri 2 (z) = ^ ( r l2 {z + Nh) R r l 2 (z) - R r l 2 (z) r 12 (z - Nh )) 

+ i (E^z + Nh) + E 1 {z- Nh) - 2E 1 (z)'"jR^(z) . 

The Fay identity in Mat(N , C)® 3 JU, [73 \8$: 

nh nh' nh' nh—h' , nh!—h nh nh nh f \ 

^ab^bc ~ f'-ac-t'ab ~r - rt f>c rt ac ) rt ab — rt a?A z a *b) ■ 


(1.48) 


(1.49) 


(1.50) 


(1.51) 


Both parts of the identity are elements of Mat (IV, C)® 3 . ft was used in | 8 ] for constructing 
higher-dimensional Lax pairs for Calogero-Moser models. Here we will prove another 
analogue of (11.1511 - in Mat (IV, C)® 2 . 

The Fay identity in Mat(N, C)® 2 : 


i&(z)i£(-t») = (1-52) 

Ncp(NH', 4-^- + S'-ft) R^ h \z+Nh') - N<p(Nh, + ft' —ft) R'^’ 1 ' (m+JVft) 


Z — 7/1 z—w y — in z—w 

+N(j){-w,^^ + h'-h)R 1 ? (w+Nh)-N<f>(-z,—— + ti-h)R 1 £ (Z + Nh '). 

The scalar analogue of this identity is obtained as follows: apply (II. 15[) (with x — h, 
y = h') to (p(h, z)4>(h', —w), and then apply (11 . 15[) once again to the obtained r.h.s.. Then 
we get the scalar analogue of r.h.s. of (11.521) . 

• Degenerated Fay identities in Mat(N, C)® 3 

RabRl = RaHab + r bc R h ac - d n R h ac , (1.53) 

+ NFl-\-z)P ac , (1.54) 

where F* h (u) — d u R' b (u ) and R bb 01 is from ill 4411 - ill .4511 

3 The identities for derivatives of i?-matrix with respect to the Planck constant and spectral parameter were 
found in [3] and [15] respectively. Authors of mm used different normalization of the i?-matrix. 
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• Degenerated Fay identities in Mat(N, C )® 2 HI. 52 1) : 

■y _ a n 

R i2(z) R 2i(~ w ) = N(f)( —— ,Nh ) (r 12 + IV/i) - 7 -i 2 (w + IVTi)) 

an — 2 ^ z — w an — y z—w 

+N$(- lr , z) RJ (z + Nh) - <o) fli 2 " (<o + iVfi) (1-55) 

2 " — 7 /) ^ — 111 
+N 2 1 ® 1 IV/i) (^(iVn) - Ei(Nh + -j^)), 

and 

< 2 (-)^ 2 i(--) = ^0^' - ^ -z) (r 12 (z + Nh) - r 12 (z + IV/i')) 

—iV0(n' - h, Nh) R h l2 h \z + Nh) + IV0(/t' - h, Nh') R h l2 n \z + Nh’) (1.56) 

+N 2 1 ® 1 <f>{h! — h , —z) (Ei(z) — Ei(z + h — h')). 

• Geometric interpretation. Due to the quasi-periodicities (ll.38|) - (ll.41l) the R- matrix have 

the following geometrical interpretation. Let V\ (V 2 ) be a rank N and degree one vector 
bundle over elliptic curve E^ with coordinate Z\ (E^ with coordinate z 2 ). Consider the 
bundle V\ K! V 2 over Er 1 '* x Let Aut-pa^N) (V\ Kl V 2 ) be the automorphism group of 
the bundle (the gauge group). The sections V(Autp GL ^ N )(Vi Kl V 2 )) depends only on the 
anti-diagonal E r of E^ x Er 2 '* with the coordinate z = Z\ — z 2 . Let E' r be the dual curve, 

h is the coordinate on E(. and V is the Poincare bundle V over S T x E(. (11.51) . Then the 

/?-matrix (11,241) is a section 

R l 2 { z ) G T ((A«tpGL(JV)(^l ^ V2)) ® V) ■ 

• Green function. Similarly to (11.191) the P-matrix can be considered as the Green function 
of (9-operator: 

BR h 12 (z) = NP 12 S 2 (z,z). (1.57) 

Properties fll.30l) - fll.48j) simply follows from their scalar counterparts except (jl.33j) which 
follows from the unitarity condition (11.271) in the classical limit (ll.30j) . Identities for derivatives 
(jl.49|) . (I1.50p were obtained in [3, Ej. Degenerated Fay identities (11.53(1 . (11.54[) in Mat(IV, C )® 3 
follows from the nondegenerated one (11.511) and local expansions (11.30(1 . (11.34]) . 

Our main interest (in this paper) is the Fay identity in Mat(IV, C )® 2 (11.521) and its de¬ 
generations (11.551) . (jl.56jl . We prove them below. The computational trick is based on the 
’’arguments symmetry” property (11.291) and the scalar Fay identities (1 1.1 51) - (1 1.1 71) . 


Painleve VI. As an application of the obtained formulae we construct higher-dimensional 
Lax pairs for the Painleve VI equation. Denote the half-periods of the elliptic curve S T as 


{D a , a = 0,1, 2, 3} = {0 , 


1 1 + r r 


} 


(1.58) 


The Painleve VI equation in the elliptic form 


is 


d 2 u 


f Jrj - 2 — U a.P ( U + • 


a =0 
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(1.59) 























































Let N be an odd (positive) integer. Consider the following pair of biock-matr icetQ: 
wm _ 1 dw / 1 <g) 1 0 \ v a ( 0 lZy 2 (u) \ 

U _ 2 d^ V 0 -! ® 1 ) ^ ivT^l V ftSiVu) 0 J 


(1.60) 


where 


and 


with 



^^“(w) = exp(2TnNhd T fl a )Ri 2 (u + iVh2 a ), 
ftSi°(-u) = exp(—2 mNhd T n a )R^(~u - NQ a ), 
TS°(m) = exp(2mNh d r Q a ) F^ 2 (u + NQ a ), 

F 2 \{~ u ) = exp(—2mNhd T il a )F 21 (—u — !Vf2 a ) 

*£(«) = $,/&(«) • 


(1.61) 


(1.62) 


(1.63) 

(1.64) 


The matrices £ Mat(2,C) <g) Mat(iV, C)® 2 . Their size equals 21V 2 x 21V 2 . The 

Painleve VI equation (11.59ft is equivalent to the monodromy preserving equation 


Tr m -{^Tn M(h) = lm ’ Mm ’ (L65) 

where the Planck constant h plays the role of the spectral parameter (see 0 ). 

For N = 1 the answer (11.601). (11.611) reDroduces the elliptic 2x2 Lax Dair DroDosed in HU. 

The Lax pair (11.601) . (II.61ft works for even IV’s as well. But the Painleve equation in this 
case has only one free constant: 


d 2 u 

dr 2 


3 



< 2=0 


( 1 . 66 ) 


2 Kronecker double series and Baxter-Belavin R- matrix 

Following idea suggested in ra we derive here the Baxter-Belavin /2-matrix as generalization 
of the Kronecker series. 

/2-matrix in Jacobi variables. Represent the elliptic curve S T (II .4[) in the Jacobi form 

C q = C/q 7, , q — e(r) = exp 2mr . 

Consider the product C q x C q with the coordinates s = e(w), t = e(z). Instead of the Kronecker 
function (J)(z,u) we consider the distribution g(s,t ) on the space of the Laurent polynomials 
C[[s -1 , £ -1 , s, £]]. For |g| < \t\ < 1 it can be represented as the series 

+n 

sM9) = Z)-^3T- (2 -1) 

n£ Z ” 

4 The coefficient l/y/—2 gives the normalization of the constants as in (11.591) . 
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If simultaneously \q\ < |s| < 1 then 

g(s, t\q) = -g + (s, t\ q) + g~(s, t\ q) , g + (s, t\ q) = ^ s*q m t n , g~(s, t\q) = ^ s l q m t n (2.2) 

z,n>0 i,n< 0 

or 

a(s,t| «) = 1 - -2— - -2-ha-(a.i) - (2.3) 

1 — 6 1 — 5 ' 

Z,71>0 

In the domain \q\ < |t| < 1 and |g| < |s| < 1 we have 


g(s,t\q)\ s= _L_ lnu ,t=^inz = (P(z,u). 

The distribution g(s,t\ q) has the properties analogous to (II .6f) ~ (II.9[i . In particular, 


(2.4) 


s(MU) =g(t,s\q). 


(2.5) 


It follows from (j2.2[i that 


g(s 1 ,t *| q) = —g(s, 1 1 q) + S(t) + 5(s) - 2 , (2.6) 

where S(s) is the distribution on the space of the Laurent polynomials 

C[t,t _1 ] = {'ip( t ) = J2i c rt l }i defined by the functional (<5,-0) = Res\ t =oi^(t) and represented by 
the formal series 

m = Y. e - < 27 ) 

nEZ 

The analog of the quasiperiodic property (II. lip is the following. The distribution g(s,t ) is a 
solution of the difference equation on t (the Green function) variable 


sg(s, tq\ q) - g(s, t\q) = S(t)-l. (2.8) 

It defines the continuation of g(s,t\ q) from the annulus |g| < |i| < 1 to C*. Due to (j2.5[i the 
similar equation can be written with respect to the s variable. 

Let 77 = e(h). The R- matrix (11.241) takes the following form in variables (s,t, rj): 

RUs)= s a * /N g{s,uj a + h)T a ®T_ a = (2.9) 

ckE Zn xZjv 

s a2/N ( ^{na 1 /N)q n ^ n+a ^ N \ n s m ) T a ® T_ a . 

aSZ^xZjv \m,n / 

It plays the role of the Green function for the difference operator 

V (A ® 1) R h 12 (sq) (A" 1 < 8 > 1 ) - R* 2 (s) = (S(s) - l)P i2 • (2.10) 


Kronecker double series 


m 


The distribution g(s, t\ q) (and (j)(z, u)) can be represented as a Kronecker double series. Consider 
the lattice in C 

W = {7 = m + nr , m, n G Z} . 
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Represent the argument u of (j>(z,u) as u = U\ + u 2 r (ui,u 2 are real), and let 

Xu ( 7 ) = e(-mu 2 +nu 1 ) 

be a character of the lattice W (xu( 7 ) : ID —> S' 1 ), parameterized by u 6 E r . The Kronecker 
double series is defined as: 

C / 1 \ Yu (7) 




* + 7 


From the definition we find that 

S(z + 1, u | r) = e(it 2 )S(z, n| t) , 

S(z + r, u| r) = e(— tti)S(z, u| r). 

It was proved in [16] that S(z,u | r) is related to the Kronecker function as 

S(z, u | r) = e(u 2 z)cp(z, u ), 

or in the Jacobi coordinates 

S(t,s\q) = i™g(s,t|g) • 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 

(2.14) 


Let us now pass to the R-matrix and describe it in terms of the Kronecker double series 
S(z,u\ T ) (|2.1ip . 

Define the lattice W by the two generators (ai/N + hi, (a 2 /N + h 2 )r), where h = hi + h 2 r, 
hi 2 G M. The corresponding character of W is 

X(m,n)(oi, h) = e ( -m{a 2 /N + h 2 ) + n(ati/N + hi)) . (2.15) 

Then the R-matrix (11,241) is dehned in terms of the Kronecker double series (12.111) as 


E X(m,n)(a, h)T a ®T_ c 
KiW = e(-h 2 z) - 

(ra,n)EZ®Z 


z + m + nr 


(2.16) 


The quasi-periodicities (ll.40p . (I1.4ip now become evident. It follows from (I2.13P that the singular 
behavior z, h —> 0 of this representation is in agreement with (ll.36p . 

We pass from Ry 2 (z) to the modihed matrix 

Ruiz) = e(h 2 z)R r l 2 (z). 

It satisfies the Yang-Baxter equation and has the quasi-periodicities 

R h i 2 {z + 1) = e{h 2 ){Q~ l ® 1 )R\ 2 (z){Q ® 1), 

RUz + t) = e(hi) (A -1 ® 1)-Ri 2 (^)(A ® 1), 

(compare with (11.401) 1. In contrast with (j 1.41 p R is not holomorphic in h and is double-periodic. 


11 






















Remark 1 The representation \2.1fy) means that the elliptic R-matrix is represented as the 
averaging of the Yang matrix z~ l P\ 2 along the lattice W twisted by the character H2.15\ ) . 

From (11.301) we also find the representation for the classical r-matrix: 


r l2 (z) = E-l(z) 1 ® 1 + Y 

m,n€(Z©Z)\( 0,0) 


E X(m,n)( a , 0) T a (g) T_ c 

OtG Z N xZjv 

z -\- m + nr 


and 


Er(z) — fcjfz) \ ^ 

mi 2 (z) = 1V ' V 1 1®1+ Y 

m,n£(Z® Z)\(0,0) 


E (2 + m + nr)x(m,n) («, 0) <g> T_ c 

ckG Zjy xZjy 

(z + m + nr) (f — r) 


3 Derivation of identities 


Proposition 3.1 The R-matrix satisfies the arguments symmetry property II 1.291) . 

Proof: Using definitions (11. 28p and (I1.23P we have 

Rfil (Nh)P 12 = JT Y TaT P ® T -u T -P TcfiNh, uj a + j-) 


a,/3 


jy 'y 1 K a,P Ta+P ® T- a -p <p a (Nh, Ul a + ^) • 


(3.1) 


a,f3 


Since K a> p = n^a+p-, the property (11.291) is equivalent to the following set of N 2 identities: 


L ^ ^ 

Jf 4, 7 ^{Nh, UJ a + ^) = <p^(z, UJ^ + h), V 76 Z 


z 

N 


x2 


(3.2) 


or 


^ _ 2, 

7 722 K U* , «( z ' u “ + fi ) = PpNft.u-i + jj), V 7 e z x2 . 


(3.3) 


The latter is verified by comparing residues. To do it we also need the relation for the sums of 
TV-th roots of 1 (it also follows from Pf 2 = 1): 




(3.4) 


Let us calculate the residue of both parts of (13.21) at h — The answer for the r.h.s. is 

obviously 5^ 1 exp(2'Kid T uj 1 z) due to (11.81) . For the l.h.s. we have: 


fi Res — Y 4, 7 Pa(Nh, u a + ^- r )= Res — Y 4 , 7 V*( Nh ~ + Jf) 


jl.l0t,QJJJ 1 ^ 2 2 / Q a \ t ATP . ^ \ 

fe s Tv exp {2-indrUJ^z) (p a {Nh, u a + —) 


(3.5) 


H3J 1 2 / 0 o x 1 B r / 9 rx x 

> J K ar/ _ u exp(27rzd T cu^)— = exp(27r*d T tu^) 
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Proposition 3.2 The R-matrix 1 1.2$ satisfies the Fay identity lil.52\) in MafiN, C)® 2 . 
Proof: Consider 


R\ 2 {z)R^ 1 (-w) = - ^ K 2 a /3 T Q+/3 ® T_ a _p tp a (z, uj a + ft) (pp(w, tup - h!) = ( 3> g\ 

a, (3 

Here we already used R^'fi—w) = —R \2 i w )- Apply the Fay identity (11.151) . then (l).•ill , and then 
(I1.15H again: 

= «a,/9 T a+/3 ® T_ a _p <p a (z - W, UJ a + ft) (f a+ p (w, U a+ p + ft - ft!) 

a,/3 

(3.7) 

- y~| V./V. + .i <8 T ,, :i ip B (w - Z,Uj J - ft') ip a+ p(z,U a+e + ft - ft') 

a,P 


= — N ^2 r 7 <g) T _ 7 ip^(Nh, uj 1 H-jy - ) ^ 7 w 7 + h — ft!) 

7 

+N Ty ® T-ry (pj(Nft!, co-y H - jj— ) (p 7 (z, a; 7 + ft — ft!) 

7 

= iv^ r 7 <g> r_ 7 ( -<t,(Nft, + ft’ - ft) ipfiw + Nft,^ + ft- ft') 

7 


(3.8) 


., , .. z — w. . Ar . z — w. 

-(pyw, ft — ft -——) </? 7 (w + A'n, a; 7 H- ——) 


N 


N 


(3.9) 


z — w 


-\-(J)(Nh , ——-h h — /i) -j - ? cj/y -I - H — h ) 


TV 


2 — 7/; ^ — 7/; 

+0(-, ft-ft 1 -—) </3 7 (~ + Nti, w 1 + —^—) 


Proposition 3.3 The R-matrices !jl.2$ and hl.Sl\ ) satisfies the degenerated Fay identities 

mm) . mm) m Mat(N, c)® 2 . 


Proof: We begin with (1 1.55 f) . Consider 

Ri 2 (z)R%i{~w) = - ^2 K l,p T a+P ® T- a -p <p a {z,u a + ft) (pp(w,up - ft). 

a,/3 


(3.10) 


Subdivide it into two parts: Yap = Y a ^~p + Y a =~p- The P ar t * s transformed as in the 
previous Proposition (via (11.1511 . then (13.311 . and then (II. 15^ again) 

^ ^ K 2 a pT a+ p ® T- a _p (p a (z ~W,UJ a + ft) tp a+ p(w , UJ a+ p) 

a^-P 

(3.11) 

- ^ K’l.pTa+fl ® T-oc-f .j lp„(w -Z,U 0 - ft) (J, U„ +/J ) 
a^-P 
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= - = -m z — T 7 g T_ 7 + Nh , a; 7 ) 

77^0 

-N*C^±,v,)'E T 7 g T_ 7 <^ 7 (u; + TV/i, (u 7 H-——) 

77^0 

(3.12) 

7 (/7 7 (z + TV A, CJ 7 ) 

77^0 

+ N 4>('^ 1 ,z)Y, T 1 ® r -7 V7(‘ + 

77^0 

By adding (and subtracting) scalar terms (lgl) to each line one obtains the first and the second 
lines of (11.551) . The input to the scalar part should be summed up together with 

^ = 1 g 1 (p a (z, UJ a + h) <Pa(-W, UJ a + h) 
a=—(3 a 

lgl ^ (p a (z — w,uj a + K) (Ei(z)—Ei(w) + Ei(h + u a ) — Ei(z — w + h + u Q )) ■ 

(3.13) 

The latter expression is transformed via (13.31) for 7 = 0 

^2,y a {z- w, u> a + h) = N(j)(Nh, —) 

a 

and its derivative (jl . 13[) . (II . 14)1 with respect to h: 

- W,CU a + k) (Ei(z - W + H + U a ) - Ei{h + (jJ a )) 

a 

= N 2 ij>(Nh, —jf~) (E 1 (Nh + ^)-E 1 (NK)\ . 

This finishes the proof of (11.551) . The identity (II .56[) can be derived similarly. Equivalently, 
(I1.56P follows from (II.55f) by using the properties (II.29p and (11.371) . ■ 


4 Higher-dimensional elliptic Lax pairs for Painleve VI 


Different types of matrix-valued Lax pairs for Painleve equations are known (see e.g. EEmm). 
In this section we construct R- matrix valued generalization of the elliptic 2x2 Lax pair suggested 
in [TTj . 


Proposition 4.1 The Painleve VI equation in the elliptic form 111.591) is equivalent to the mon- 
odromy preserving equation l\1.65\) with the Lax pair lil.60\) - \1.64\) and the elliptic R-matrix 
jXIg] ) for odd N. 


Proof is similar to the one given in S3 for the scalar (TV = 1) case. First, notice that fj^L(K) = 
2 'r-d u L(H ) + d T L(h), where the last term is the derivative by explicit dependence on r. It is 
canceled out by ^^M(^) due to the heat equation (11.481) 2m d T lZ r ]fff{u) = d^J-^iu). 
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Denote 


L a = 


n 


n h 2l \-u) 


h,a, 

12 1 


U 


M a = 


0 

21 (~~ U ) 


r 




(4,1) 


The main statement which we need to verify is that for a ^ b 

[L a , M b ] + [L b , M a ] = 0 , (4.2) 

i.e. the input to comes only from [L a ,M a ]. Indeed, it follows from the unitarity 

condition (11,271) that 

Ui 2 (u)n b [ a (-u) = R h l2 {u + NQ a )R^(-u - NQ a ) = N 2 (p(Nh) - p(u + NQ a )). (4.3) 

Differentiating (14.31) with respect to u we get 

^ 2 » 7 Z h 2 i(-u) - = —N 2 p'(u + NQ a ). (4.4) 


This identity provides the equation of motion. Notice that in order to have all four constants N 
should be odd since p'(u + NQ a ) = p'{u-\-VL a ) in this case. If N is even then p'(u + NQ a ) = p'(u), 
and we have only one constant as in (11.661) . 

To prove (j4.2[) let us recall that in the scalar case this followed from 


U + Q a )f b (h, -u - Q b ) - f b (h, u + Q b )(p a (h, -u - Q a ) 
tp b (h, U + Q b )f a (h, -u - Vt a ) - f a (h, U + VL a )p b {h, -U-Vt b ) = 0 , 


(4.5) 


where 

f a (z,u + Q a ) = exp(2v Tid r tt a h)d w (j)(h,w) | w=u+no 

is the scalar analogue of T\ 2 (w) • The identity (14.51) appears from (11.161) and (ll.lOI) - (ll.lip as 
follows: 

p a (h, u + n a )p b (h, u - n b ) + p b (h, u + n b )p a (h, -u - n a ) = 

(4.6) 

^a+b{^i ^a + ^ b ) {2Ei(h) — Ei(h + Q a — Q b ) — Ei(h + Q b — D a )). 

The r.h.s. of (j4.6[) is independent of u. The derivative of (|4.6[) with respect to u gives (14.51) . 
Similarly to (14. 6 j) it follows from the degenerated Fay identity (1 1.55 j) that 

7 z h £{u)n h *{-u) + n h *{u)n h £{-u) 

(4.7) 

= iv 2 1 <g> i ip a+b (Nh, n a + n b ) (2 e 1 (nh) - e^nh + n a - n b ) - e^nh + n b ~ n a )). 

It can be verified directly using (ll.lOI) - (ll.lll) which can be re-written as 


( t>(z, w + D a ) = exp(— 2mzd T Q a )(p(z, w — D a ). 

The r.h.s. of (14.71) is scalar and independent of u. The derivative of (14.71) with respect to u gives 
E^ a (u)n H 2 f(-u) - R h 12 a (u)E b i b (-u) + E h l2 b {u)H h 2 ?{-u) - n% b (u)E b { a (-u) = 0 . (4.8) 

This identity underlies (14.21) . ■ 
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